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Transport theory for a two-flavor color superconductor 



o 
o 

(N 
C 

ON 



(N 
> 

On 
O 

m 
o 



43 

Oh! 
i 

Oh! 
<D 
43 



X 



Daniel F. Litim and Cristina Manuel 
Theory Division, CERN, CH-1211 Geneva 23, Switzerland. 

QCD with two light quark flavors at high baryonic density and low temperature is a color super- 
conductor. The diquark condensate partially breaks the SU(3) gauge symmetry down to an SU(2) 
subgroup. We study thermal fluctuations of the superconductor for temperatures below the gap A. 
These are described by a simple transport equation, linked to a quasiparticle behavior of the thermal 
excitations of the condensate. When solved in the collisionless limit and close to equilibrium, it gives 
rise to the "hard superconducting loop" (HSL) effective theory for the unbroken SU (2) gauge fields 
with momenta t< A. This theory describes Debye screening and Landau damping of the gauge 
fields in the presence of the diquark condensate. We also explain how our effective theory follows 
to one-loop order from quantum field theory. Our approach provides a convenient starting point for 
the computation of transport coefficients of the two-flavor color superconductor. 

PACS numbers: 12.38.Mh, 24.85. +p 



Quantum chromodynamics (QCD) under extreme con- 
ditions displays a very rich phase structure. At low 
temperature and high baryonic density, quarks form 
Cooper pairs due to the existence of attractive interac- 
tions amongst them. The diquark condensate modifies 
the ground state of QCD and leads to the phenomenon 
of color superconductivity This phase is typically 

characterized by the Anderson-Higgs mechanism and an 
energy gap associated to the fermionic quasiparticles. It 
is expected that this state of matter is realized in com- 
pact stars. Studying the properties of cold dense quark 
matter may explain some basic features of neutron stars 
and could even lead to the prediction of new stellar ob- 
jects (see for a recent review on the subject). 

A lot of progress has been made in the microscopic 
understanding of color superconductivity. Important 
open questions concern macroscopic properties of dense 
quark matter. Transport theory is known to provide 
a very efficient framework for the study of low-energy 
and long-range fluctuations in a medium. In this Letter, 
we present a transport equation for the gapped quasi- 
particle excitations of the two-flavor color superconduc- 
tor. This is a convenient starting point for both the con- 
struction of low energy effective theories and the compu- 
tation of macroscopic observables within the kinetic or 
hydrodynamical limits. 

The global properties of a color superconductor de- 
pend dramatically on the number of quark flavors Nf 
participating in the condensation. Here we shall discuss 
the case of two flavors Nf = 2, which exhibits a very 
rich structure. The diquark condensate breaks the color 
SU(3) group into an SU(2) subgroup. As a result, five 
gluons acquire masses through the Anderson-Higgs mech- 
anism, while three of them remain massless. Also not all 
the quarks attain a gap. In the high-density limit, and 
because of asymptotic freedom, the gap and the gluon 
masses can be computed reliably from first principles 
]^-|Tl|]. A well-defined hierarchy of scales then appears: 
at weak coupling g <C 1, the gap A is exponentially 



smaller than the gluon masses ~ g/i and the chemical 
potential fj,. 

The low energy physics of a two-flavor color super- 
conductor is dominated by its light degrees of freedom. 
At vanishing temperature, these are the massless gauge 
bosons, the gapless quarks and a (pseudo-) Goldstone 
boson, similar to the r] meson. However, the gapless 
quarks and the r\ meson are neutral with respect to the 
unbroken SU(2) subgroup. In turn, the condensate, al- 
though neutral with respect to the unbroken SU(2), po- 
larizes the medium since their constituents carry SU(2) 
charges. Hence, the dynamics of the light SU (2) gauge 
fields differs from the vacuum theory. This picture has 
recently been introduced by Rischke, Son and Stephanov 
12 1 . Their infrared effective theory for momenta k <C A 
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where Ef = Fffc and Bf = \e ijk F9 k are the SU(2) 
electric and magnetic fields. The constants e and A 
are the dielectric susceptibility and magnetic perme- 
ability of the medium. To leading order, A = 1 and 
e = 1 + g 2 /x 2 /(187T 2 A 2 ) Q. As a consequence, the ve- 
locity of the SU (2) gluons is smaller than in vacuum. 
This theory is confining, but the scale of confinement is 
highly reduced with respect to the one in vacuum with 



A! 



Q( ~ Aexp (-M*ZJ±) Q. Due to asymptotic free- 
dom, it is expected that perturbative computations are 
reliable for energy scales larger than Aq CD . 

At non-vanishing temperature, thermal excitations 
modify the low energy physics. The condensate melts at 
the critical temperature T c ss 0.567Ao @] (Ao is the gap 
at vanishing temperature). We restrict the discussion to 
temperatures within Aq CD <C T < T c , which provides 
the basis for the perturbative computations below. In 
this regime, the main contribution to the long distance 
properties of the SU(2) fields stems from the thermal ex- 
citations of the constituents of the diquark condensate. 
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The thermal excitations of the massless gauge fields con- 
tribute only at the order g 2 T 2 and are subleading for 
sufficiently large [i. Those of the gapless quarks and of 
the 77 meson do not couple to the ST/ (2) gauge fields. 

The thermal excitations due to the constituents of 
the diquark condensate display a quasiparticle structure. 
This implies that they can be cast into a transport equa- 
tion. To that end, and working in natural units ks — h — 
c = 1, we introduce the on-shell one-particle phase space 
density f(x,p,Q), x M = (t,x), describing the quasipar- 
ticles. The distribution function depends on time, the 
phase space variables position x, momentum p, and on 
SU(2) color charges Q a , with the color index a = 1,2 
and 3. The quasiparticles carry SU(2) color charges sim- 
ply because the constituents of the condensate do. The 
on-shell condition for massless quarks m q — relates 
the energy of the quasiparticle excitation to the chemical 
potential and the gap as pa = e p = \J (p — p) 2 + A 2 (T). 
The gap is both temperature and momentum-dependent. 
From now on, we can neglect its momentum dependence 
which is a subleading effect. The velocity of the quasi- 
particles is given by 



de p 
dp 



b — mI 



V(p-/x) 2 + A2(T) 



(2) 



Transport equations similar to Eq. (0) have been 
known previously: (i) The analogue of Eq. (||) at high 
temperature for (j, = A = (and SU(N) color charges 
Q a ) has been introduced in [[13[ for the description of the 
quark-gluon plasma at high temperature. In the Vlasov 
approximation and for C[f] = 0, it reproduces the hard 
thermal loop (HTL) effective theory |l4j to leading order 
in g [^5). (ii) A formalism which allows to go systemati- 
cally beyond the Vlasov approximation has been derived 
as well |l6|]. This includes the derivation of the collision 
term C[f]. (iii) At high density and T = A = 0, the cor- 
responding transport equation leads to the hard dense 
loop (HDL) effective theory for C[f] = 0, and to leading 
order in g pJ^ . (iv) In the non-relativistic limit for the 
Abelian case, our transport equation reduces to the one 
for a BCS superconductor 

In the remaining part of the Letter, we study the 
collisionless dynamics C[f] = of the color supercon- 
ductor close to thermal equilibrium and to leading or- 
der in the gauge coupling. Consider the distribution 
function f(x, p,Q) = f cq (p ) + gf (1) (x,p,Q). Here, 
f eq '(Po) = l/( ex P( e p/T 1 ) + 1) is the fermionic equilib- 
rium distribution function and gf^(x,p,Q) describes a 
slight deviation from equilibrium. For convenience, we 
also introduce the color density 



and depends on both the chemical potential and the gap. 
For A = 0, the quasiparticles would travel at the speed 
of light. However, in the presence of the gap A ^ 0, their 
propagation is suppressed, v p = |v p | < 1. 

The one-particle distribution function f(x, p, Q) obeys 
a very simple transport equation, given by 



Dt + v p --D- gQ a (E a + v, 



op 



f = C[f}. (3) 



Here, we have introduced the short-hand notation D^f = 
[dfj, — ge abc Q c A^d^]f for the covariant derivative acting 
on /. The first two terms on the l.h.s. of Eq. <M) combine 



to a covariant drift term v^D^ 7 where vp 1 = (l,v p ) and 
= (D t , D). The terms proportional to the color elec- 
tric and magnetic fields provide a force term. The r.h.s. of 
Eq. (ffl contains a (yet unspecified) collision term C[f]. 

The thermal quasiparticles carry an SU (2) charge, and 
hence provide an SU(2) color current. It is given by 



■TO = 9 



E / jffisdQv$Q*f(x,p,Q) . (4) 



Below, we simply omit a species or helicity index on /, 
as well as the explicit sum over them. The color mea- 
sure is normalised J dQ — 1, and obeys J dQQ a = 
and J dQQaQb = C^ah, where C2 denotes the quadratic 
Casimir (C2 = \ for quarks in the fundamental). The 
color current Eq. (||) is covariantly conserved for C[f] — 
0. For C[f] / a covariantly conserved current implies 
certain restrictions in the form of the collision term. 



J a (x, p) = g dQQ a f(x, p, Q) , 



(5) 



from which the induced color current of the medium 
Eq. (g) follows as JP{x) = J j^vPJ a {x,p). Expand- 
ing the transport equation ([|) to leading order in g, and 
taking the two helicities per quasiparticle into account, 
we find the transport equation for the color density as 



[D t + v p • D] J(x, p) = g 2 N f v p • E(z) 



df e 



dc. 



The solution of the transport equation reads 
JP(x) = g 2 N f f^y v ^ GabVp ■ Kv:,/:'''" 



(2tt) 3 



dtr. 



(G) 



(7) 



with the Greens function G a b = (y\^/( v p • D)\x) a b- Af- 
ter having solved the transport equation, the relevant 
information concerning the low energy effective theory is 
contained in the functional J[A]. Notice that the above 
derivation is analogous to the derivation of the HTL 
and HDL effective theories from kinetic theory Jl4],[l5],[l7) . 
Owing to this resemblance, we call the diagrams which 
are derived from Eq. (Q) as hard superconducting loops 
(HSL). The HSL effective action follows from Eq. (0) by 
solving J[A] = -6T HSh [A]/5A for Fhsl[A], and all HSL 
diagrams can be derived by performing functional deriva- 
tives to the effective action (or the induced current). We 
thus reach to the conclusion that the low energy effec- 
tive theory for a two-flavor color superconductor at finite 
temperature reads S^ s = S^f + Thsl to leading order 
in g. This theory is effective for modes with k -C A. 
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Let us have a closer look into the induced cur- 
rent, which we formally expand as J£[A] = W^ V A^ + 
^rffi A%AP + ... in powers of the gauge fields. The most 
relevant information on the thermal effects is contained 
in the thermal polarization tensor II °* . Using Eq. (Q), 
we find 



cq 



d 3 p df 
(2tt) 3 ~d7 p 



g »o g »o_ kQ ^L 



It obeys the Ward identity fc^IT^fc) = 0. With retarded 
boundary conditions fco ~ * fco + *0 + , the polarization ten- 
sor has an imaginary part, 



ImlC(fc) = 6 ab irg 2 N f k 



d 3 p df oc 
(2tt) 3 ~dT p 



which corresponds to Landau damping. Performing the 
angular integration, we obtain for the longitudinal and 
transverse projections of the polarization tensor 



n L (fc) = 



n T (fc) = 



9 2 N f 
2tt 2 

In 



°, 2 df e 
dpp -T- 
o de } 

fco 



, 1 fco 
1 - -- — X 

2 kv n 



kv 



kv 



fco 

g^kj 
4tt 2 k 2 

fco + kv p 



°, 2 #° q 
d PP -J— 

de n 



iirQ{k 2 vl - fcg) 
kv. 



(8a) 



In 



fco — kv p 



2„,2 



fco 

kv, 



Mr6(fcV - ko) 



(8b) 



where is the step function. We first consider the real 
part of the polarization tensor. From Eq. (g), and in 
the limit fco ~* 0, we infer that the longitudinal gauge 
bosons acquire a thermal mass, the Debye mass, while 
the transverse ones remain massless. The (square of the) 
Debye mass is given by 



' D 



9 2 N f 
2tt 2 



dpp 



2 df cq - 
de n 



M 2 In I = - 



A T 



(9) 



For convenience, we have factored-out the Debye mass 
M of the ultradegenerate plasma in the normal phase, 
M 2 = g 2 Nf[i 2 /2ir 2 . The dimensionless functions 



(- 








l) 


P 2 l 



• df 



cq. 



dCr, Vp 



(10) 



obey > > for all n due to v p < 1. Equality 

holds for vanishing gap. For the physically relevant range 
of parameters T < A <C [i, the functions /„ are <C 1. In 
particular, it is easy to see that Z„(oo,0) = 0: there is 
no Debye screening for the SU (2) gluons at T = in the 
superconducting phase. In the limit where A/T 3> 1, and 
to leading order in T//i 1, the Debye mass reduces to 



\/ 2 A/27r^exp(-A/T). 



(11) 



The dispersion relations for the longitudinal and trans- 
verse gluons follow from the poles of the corresponding 
propagators, 

efc 2 - ReLI L (fc , fc)U =a, i( fc) = , (12a) 

e fc 2 -ifc 2 + Ren T (fc 0) fc)U 0=a , T(fe) =0 , (12b) 

where e and A have been introduced in Eq. (|l]). The 
plasma frequency uj p i follows from Eq. (12) as 



1 



= —M'h -,- 



3e 



A T 



T /i 



(13) 



For generic external momenta the dispersion relations 
can only be solved numerically. In turn, if the spatial 
momenta are much smaller than the plasma frequency 



ers of k 2 /u 2 y as 



fc <C lu p \. solutions to Eq. (12) can be expanded in pow- 



u,|(fc) 



Cjf^fc) 



-'pi 



"pi 



1 + 



3/4 fc 2 

5 72 ^ 



k 4 



pi 



eA 



:.2 



1 lh\ fc 



-pi 



"pi 



(14a) 
(14b) 



Let us now consider the imaginary part of Eq. (g) , which 
describes Landau damping. Since v p < 1, we conclude 
that Landau damping only occurs for fc 2 < fc 2 . Hence, 
plasmon and transverse gluon excitations are stable as 
long as LOL,T(k) > fc. Furthermore, we notice that the 
imaginary part of Eq. (|8|) is logarithmically divergent: 
the quasiparticle velocity vanishes for momenta close to 
the Fermi surface, which is an immediate consequence of 
the presence of a gap, cf. Eq. (|^). This divergence does 
not appear in the real part, because the logarithm acts 
as a regulator for the l/v p factor. To leading order in 
T//i, and in the region of small frequencies fc 2 -C fc , we 
find at logarithmic accuracy, and for all values of A/T, 



ImII L = -2irM 



fco A 



2 « 



ln(fc/fc ) 



ImlLr = ttM 2 ^ 



;fco 
fc 



fc T(e A / T + l)(e- A / T + l) ' 
1 



fc 2 A 



e A / T + 1 

In fc/fco 



J pT(e A / T + l)(e- A / T + l) 



(15a) 



(15b) 



For small frequencies, Landau damping is dominated by 
the logarithmic terms, which are proportional to the gap. 
Once the gap vanishes, subleading terms in not dis- 



played in Eq. (15), take over and reduce Imll to known 
expressions for the normal phase. 

Finally, we explain how the polarization tensor, as ob- 
tained within the present transport theory, matches the 
computation of IP 1 ' for external momenta fco,fc <C A 
to one-loop order from quantum field theory. The one- 
loop gluon self-energy for a two-flavor color supercon- 
ductor has been computed by Rischke, and the po- 
larization tensor for the unbroken SU(2) subgroup is 
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given in Eq. (99) of pM. It contains contributions from 
particle-particle, particle-antiparticle and antiparticle- 
antiparticle excitations. The particle-antiparticle con- 
tribution to I! 00 and n 01 at low external momenta, 
and the antiparticlc-antiparticle excitations are sublead- 
ing. The particle-particle contributions divide into two 
types. The first ones have poles for gluonic frequencies 
ka = ±(e p + e p _fc) and an imaginary part once ko ex- 
ceeds the Cooper pair binding energy 2A. These terms 
are related to the formation or breaking of a Cooper 
pair, and suppressed for low external gluon momenta. 
The second type of terms, only non-vanishing for T ^ 0, 
have poles at fco = ± (e p — e p -fc). For I- « A we ap- 
proximate it by fco ~ ^TJJT ' k- The prefactor, a differ- 
ence of thermal distribution functions, is approximated 
by / cq (e P ) - / cq (e P -fc) « ^ • k*g. After simple al- 
gebraic manipulations we finally end up with the result 
given above. We conclude that this part of the one-loop 
polarization tensor describes the collisionless dynamics of 
thermal quasiparticles for a two-flavor color superconduc- 
tor. The same type of approximations can be carried out 
for n y to one-loop order. There, apart from the HSL con- 
tributions, additional terms arise due to particle-particle 
and particle-antiparticle excitations, cf. Eq. (112) of |Tl]| . 
We have not evaluated these terms explicitly. However, 
we expect them to be subleading or vanishing, as oth- 
erwise the Ward identity fc M n^(fc) = is violated. For 
T = 0, this has been confirmed in jljj ]. 

Summarizing, we have provided a transport equation 
for the gapped quarks of two-flavor color superconduc- 
tors. Its simple structure is based on the quasiparticle 
behavior of the thermal excitations of the condensate, 
in consistency with the underlying quantum field the- 
ory. We have constructed a low temperature infrared 
effective theory of the superconductor. To leading or- 
der, we found Landau damping, and Debye screening of 
the chromo-electric fields. Beyond leading order, chromo- 
magnetic fields are damped because they scatter with the 
quasiparticles. The damping rate is related to the color 
conductivity. It should be possible to compute the rate 
from the transport equation ([}]), amended by the relevant 
collision term. The latter can be derived, for example, 
using the methods developped in Jig] , 

We have not discussed the transport equations for gap- 
less quarks or the rj meson, because they do not carry 
SU(2) charges. However, their excitations are light com- 
pared to the gapped quasiparticles, and dominant for 
other transport properties such as thermal and electri- 
cal conductivities or shear viscosity. The corresponding 
set of transport equations will be discussed elsewhere. 

It would be very interesting to study the transport 
equations in a three-flavor color superconductor |L9) . For 
Nf = 3 the quark-quark condensate breaks the SU(3) 
gauge group completely, as well as some global flavor 
symmetries. Transport phenomena should then be dom- 
inated by the Goldstone modes associated to the breaking 
of the global symmetries. The corresponding transport 



equations will be substantially different for the two and 
three flavor case. 
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